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Discovery of 7-ray emission from the Broad Line Radio 
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ABSTRACT 

We report the discovery of high-energy 7-ray emission from the Broad Line Radio 
Galaxy (BLRG) Pictor A with a significance of ~ 5.8c (TS= 33.4), based on three 
years of observations with the Fermi Large Area Telescope (LAT) detector. The three- 
year averaged E > 0.2 GeV 7-ray spectrum is adequately described by a power- 
law, with a photon index, T, of 2.93 ± 0.03 and a resultant integrated flux of F 7 = 
(5.8 ± 0.7) x 10~ 9 ph cm^s" 1 . 

A temporal investigation of the observed 7-ray flux, which binned the flux into 
year long intervals, reveals that the flux in the third year was 50% higher than the 
three-year average flux. This observation, coupled with the fact that this source was 
not detected in the first two years of Fermi-L AT observations, suggests variability on 
timescales of a year or less. 

Synchrotron Self-Compton modelling of the spectral energy distribution of a 
prominent hot-spot in Pictor A's western radio lobe is performed. It is found that 
the models in which the 7-ray emission originates within the lobes, predicts an X- 
ray flux larger than that observed. Given that the X-ray emission in the radio lobe 
hot-spots has been resolved with the current suite of X-ray detectors, we suggest that 
the 7-ray emission from Pictor A originates from within its jet, which is in agreement 
with other 7-ray loud BLRGs. This suggestion is consistent with the evidence that 
the 7-ray flux is variable on timescales of a year or less. 

Key words: galaxies: active - galaxies: Broad Line Radio Galaxies - galaxies: individ- 
ual (Pictor A = PKS 0518—458) - galaxies: jets - gamma rays - radiation mechanisms: 
non thermal 



1 INTRODUCTION 

The unified model of Active Galactic Nuclei (AGN) at- 
tributes the bright nucleus at the center of some galaxies to 
a super-massive black hole-accretion disk system surrounded 
by regions of broad line and narrow line gas clouds. Within 
this unified model, the difference between radio-loud (RL) 
and radio-quiet (RQ) AGN is explained by RL AGN possess- 
ing a relativist ic jet emanating a p proximately perpendicul ar 
from the disk (jAntonuccil l|l99of ): lurrv fc Padovanl l|l995f )h 
The early success of the unified AGN model was its abil- 
ity to explain the observational characteristics of the differ- 
ent RL AGN sub-classes primarily by the size of the angle 
between the relativistic jet and the observer's line of sight. 
However, while the unified model works well to first order, 
there is growing evidenc e that it is an oversim plification (e.g. 
lAharonian et all l|2007l ): iKharb et all (|2010l )) and there are 
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still several key fundamental questions that remain to be 
answered. 



One of these questions is why relativistic jets 
only appear to be present in 10-20% of AGN. This 
question is often included in the broader issue of 
how accretion and ejecta are linked in AG N (e.; 
Blandford fc Znaiekl dl977T); iBlandford fc Pavnel <fl982 



Sikora. Stawarz fc Lasotal l|2007t ): lMcNamara et all lj201lh h 
To address this issue, the accretion disk-relativistic jet con- 
nection needs to be studied across a wide variety of RL 
AGN types. Broad Line Radio Galaxies (BLRGs) appear 
to be vital in such a study. With the line of sight orien- 
tated at an intermediate angle relative to the relativistic 
jet, BLRGs exhibit both jet-related and disk-related photon 
emission processes. BLRGs, therefore allow us to investigate 
how the processes of accretion (disk) and ejection (jet) relate 
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to each other, and as such, represent a powerful diagnostic 
tool through which to understand this disk-jet connection^. 

The successful launch of the Fermi 7-ray Space Tele- 
scope affords us an ideal opportunity to investigate the 
accretion disk-jet connection in AGN. The ability of the 
Fermi-liAT detector to scan the entire 7-ray sky ev- 
ery three hours allows us to search for 7-ray emis- 
sion from AGN unbiased by activity state or AGN sub- 
class. Recently, the two year LAT AGN catalogue has 
been released (2LAC), w ith the 2LA C clean sample con- 
sisting of 885 sources (|Abdo et all (|2011al )). Of these 
sources, the vast majority are of the blazar sub-class, 
however there is a growing number of other AGN sub- 
classes, including radio galaxies and Narrow Line Seyfer t 
1 galaxies, de tected as 7-ray source s (lAbdo et al.l J2009al): 
lAbdo et all j2009cD; lAbdo et all (l2009dh; lAbdo et all 



j2009eD; Kataoka et all J201C|); (Foschini et all (|2011al l; 
iFoschinil l|2011bl ): iBrown fc Adams! (|201lh ). This increase 
in the number of different AGN sub-classes detected as 7- 
ray sources is extremely important for the phenomenolog- 
ical studies required to investigate the connection between 
accretion and ejection. 

To date, two BLRGs h ave been d e tected as 7-ray 
sources, 3C 120 and 3C 111 (|Abdo et all (|2010bl )). To in- 
vestigate the_£Ojmection betw een accretion and ejection in 
BLRGs, IKataoka et all (|201lh were recently able to suc- 
cessfully model the broad-band spectral energy distribu- 
tion (SED) of both 3C 111 and 3C 120 with a combined 
model incorporat i ng ac cretion disk emission, taken from 
iKoratkar fc Blaea (1 1999T) . and relativistic jet emission, taken 
from Soldi et al. ( 20081 '). Kataoka et al. also applied this com- 
bined disk-jet model to the two-year data set of a sample of 
18 X-ray bright BLRGs. While this study did not reveal 
any new 7-ray emitting BLRGs, it did find that five of the 
BLRG sample have flux upper limits that we re very close to 
the flu x level predicted by the disk-jet model (IKataoka et al.l 
|201ll )ffl Motivated by the proximity of the upper limits to 
the model prediction, we utlised a larger three-year Fermi 
data set to search for 7-ray emission from the most promis- 
ing BLRG candidates. In total four BLRGs sources were 
studied. While three of these did not show significant 7-ray 
emission, our analysis revealed that Pictor A is a high-energy 
7-ray source. 

Pictor A (z= 0.035, lEracleous fc Halpernl |2004l ); 
iLiu fc Zhangl (|2002r T) is an archetyp i cal Fa naroff-Riley II 
radio galaxy (FRII. iFanaroff fc Rilevl (11974 )), with the to- 
tal radio emission dominated by the jet termination re - 
gion observed as lar ge radio lobes (e.g. iPerlev et all (| 1997T ) : 
iTingav et al.l (|2000|)). From optical spe ctra, Pictor A is de- 
fined as a BLRG ( Danziger et al.l (|l977T l). Interestingly Pic- 
tor A's host galaxy has recently been fo und to possess a 
disk-like morphology (jlnskip et al.l (|2010l )). indicating that 
Pictor A is hosted by a spiral galaxy, and thus is one of 
the rare examples against the AGN jet-ellipical host galaxy 



1 It should be noted that while the 7-ray emission processes 
of radio galaxy jets are thought to be somewhat different to 
those of blazars, BLRGs still allow us to investigate how the 
global jet emission is related to t he accretion disk pro cesses 
dGeorganopouIos fc Kazanasl J2003h ; iGhisellini et al.l 1 120051) ). 

2 One of these 5 BLRGs was Pictor A, which hinted at 7-ray 
emission with a TS value of 20 from two years of observations. 



paradigm fsee lFoschiml (|2011bl ) for more details). Due to its 
close proximity, Pictor A affords us an excellent opportu- 
nity to study not only the disk-jet connection but also the 
termination shock environment of the radio lobefl It is not 
surprising then that Pictor A has been extensively observed 
at many wavelengths. 

Radio observations have found Pictor A to pos- 
sess mildly relat ivistic jets containi ng subluminal mov- 
ing components (|Tingav et al l (|2000l )). Very Long Base- 
line Array (VLBA) observations have also shown com- 
plex radio structure within Picto r A's radio l obes in 
the form of compact hot-spots |Tingav et alj ((2008)). 
These remarkable hot-spots, along with the relativistic 
jet, have been studied with other observatories, including 
Spitzer, Hubble Space Telescope, XMM and Chandra (e.g. 
iThomson et all (Il995l); ISimkin et all (Il999l); IWilson et all 
j200ll);lHardcastle fc Crostonl d2005l) ; [ Grandi et all J2006l); 
iMiglipri et all l|2007t ); ITingav et all IMPOST ); iMalkan et all 
( 20 111 )). These observations have found structured magnetic 
fields within both the relativistic jet and radio lobe hot- 
spots. Synchrotron self-Compton (SSC) modelling of these 
broad-band observations predicts the possibility of 7-ray 
emission through inverse-Compton (IC) scattering of the 
synchrotron photons on the relativistic electrons. However 
the predicted flux level is somewhat lower than the cur- 
rent 7-ray upp er limits from bot h the Fermi and HES S 
collaborations (I Zhang et al.l (|2009h : IKataoka et~aH l|201ll) ; 
lAharonian et all (|2008f) ). 

While Zhang et al. predicted 7-ray emission origi- 
nating from the conditions in the radio lobe hot-spots, 
the high-energy emission from BLRGs is often attributed 
to the relativistic jet v ie wed at interm e diate angles 
(I Grandi fc Palumbol ((20071 ); ISambruna et all (|2009l )). The 
question then arises as to whether the Pictor A 7-ray emis- 
sion, that we have discovered, is associated with the radio 
lobe hot-spots or with the relativistic jet. To address this 
question we perform SED modelling of a prominent hot- 
spot in the western radio lobe of Pictor A using archival 
data at other photon energies. We take two approaches. In 
our first approach we tune the various parameters available 
in the SSC model so as to best fit the SED of the radio 
to optical data and the 7-ray data. We find that it is not 
possible to simultaneously fit these photon energy bands as 
well as the X-ray data with the X-ray flux over-predicted 
compared to observations. We then take the alternative ap- 
proach of optimising the fit of the radio to optical data and 
the X-ray data. We find a good fit to this data but the pre- 
dicted flux at 7-ray energies is significantly below the flux 
that we have observed from Pictor A. This suggests that if 
the 7-ray emission is SSC in origin, it is not produced in 
the radio lobe hot-spot region where the other wavelength 
photons are observed to originate. 

The paper outline is as follows: In §2 we describe the 
Fermi-LAT observations and data analysis routines used in 
this study. Pictor A's 7-ray characteristics are given in §3. In 
§4 we present our broad-band spectral energy distribution 



3 Assuming a ACDM cosmology and Ho = 71 km s 1 Mpc 1 
llKomatsu et alj l|2009l) ). Pictor A has a luminosity distance of 
about 148 Mpc. At this distance 1 arcsecond corresponds to ap- 
proximately 700 parsecs. 
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modelling using archival data at other photon energies, with 
the discussions on the spectral energy distribution modelling 
given in §5. The conclusions are given in §6. Throughout 
this paper, a A cold dark matter (ACDM) cosmology was 
adopted, with a Hubble constant of Ho = 71 km s _1 Mpc -1 , 
f^m = 0.27 and J1a = 0.73 as der ived from Wilkinson M i- 
crowave Anisotropy Probe results (|Komatsu et al.l |200^)). 



2 FERMI-LAT OBSERVATIONS AND DATA 
ANALYSIS 

The LAT detector aboard Fermi, described in detail by 
lAtwood et ail (|2009l V is a pair-conversion telescope, sensi- 
tive to a photon energy range from below 20 MeV to above 
300 GeV. With a large field of view, ~ 2.4 sr, improved an- 
gular resolution, ~ 0.8° at 1 GeV9, and large effective area, 
~ 8000 cm 2 on axis at 10 GeV, Fermi-LAT provides an or- 
der of magnitude improvement in performance compared to 
its EGRET predecessor. 

Since 2008 August 4, the vast majority of data taken by 
Fermi has been performed in all-sky-survey mode, whereby 
the Fermi-LAT detector points away from the Earth and 
rocks north and south of its orbital plane. This rocking mo- 
tion, coupled with Fermi-LAT's large effective area, allows 
Fermi to scan the entire 7-ray sky every t wo orbits, or ap- 
proximately every three hours (|Ritj (|2007r i). 

The data used in this analysis comprises of all all-sky- 
survey observations taken during the first ~three years of 
Fermi operation, between 2008 August 4 and 2011 Au- 
gust 8, equating to a mission elapse time (MET) interval 
of 239557417 to 334486765. In accordance with the PASS 
7 criteria, a zenith cut of 100° along with a rock-angle cut 
of 52° was applied to the data to remove any cosmic ray 
induced secondary 7-rays from the Earth's atmosphere. All 
'source' class eventqj within a 6° radius of interest (ROI) 
were considered in the 0.2 < F 7 < 300 GeV energy range, 
with the lower energy cut being conservative to reduce the 
systematic errors. 

Throughout this analysis, the most recent Fermi-LAT 
Science Tools, version v9r23pl, were used in conjunction 
with instrument response functions P7SOURCE_V6. Fur- 
thermore, the most recent galactic, gal_2yearp7v6_v0.fits, 
and extragalactic, iso_p7v6source.txt, diffuse models were 
utilised during spectral fitting and upper limit calculations. 



3 7-RAY CHARACTERISTICS 

To investigate the 7-ray properties of Pictor A, we utilized 
the unbinned maximum likelihood estimator of the Fermi- 
LAT Science Tools' GTLIKE routine. GTLIKE allows us to 

4 Below 10 GeV photon energy, the 68% containment an- 
gle of the photon direction is approximately given by 8 ~ 
0.8°(.E 7 /GeV) -0 - 8 , with the 95% containment angle being less 
than 1.6 times the angle for 68% containment. 

5 'Source' class events equates to EVENT_CLASS=2 in the PASS 
7 data set. Source class events optimise the 7-ray detection effi- 
ciency as a function of residual charged-particle contamination 
and is best suited for localised sources as opposed to transient or 
diffuse sources. 



fit the data with a series of both point and diffuse sources of 
7-rays. The model used to calculate the likelihood of 7-ray 
emission from Pictor A was a combination of the most recent 
galactic and extragalactic diffuse models, as given in §2, and 
all point sources within an 8° radius region of interest (ROI) 
centered on Pictor A. Each point source was modelled with 
a power- law spectrum of the form dN/dE — AxE~ F with 
both parameters, photon index, P, and normalisation, A, 
free to vary. The co-ordinates for each point source within 
the RO I were taken from the 2nd Fermi-LAT Catalogue 
(2FGL; IXbdo et all l|2uTlal )L In addition, a 10° radius sky 
map of the three year data set was used to search for any 
other sources in the ROI that were not present in the 2FGL. 

As a further check, the spatial distribution of the test 
statistic[f| values obtained from the GTLIKE fit was also 
utilised to identify any other possible sources in the ROI 
that were not present in the 2FGL. The TS distribution in- 
dicated the presence of two TS peaks within 1° of Pictor 
A. One of these peaks has a TS value of 76.6 and is spa- 
tially coincident with the blazar BZQ J0515— 4556. Indeed, 
BZQ J0515— 4556 is present in the 2FGL and therefore was 
already included in the model file. 

The other TS peak is separated from Pictor A by 
~ 0.9°, located at the position (QJ2000, 5j20oo = 81.3°, 
-45.9°). With a TS value of 13.8, or ~ 3.7cr, this TS peak 
suggests the presence of another 7-ray source which is not in 
either the 1FGL or 2FGL catalogues. However given that the 
TS value is below the detection threshold of 25, this emis- 
sion is not yet significant enough to be formally classified as 
a 7-ray source. 

Another Fermi-LAT Science Tool, GTFINDSRC, was 
applied to the 8° ROI around Pictor A to localise the origin 
of the 7-ray emission. Using the same combined diffuse and 
point source model that was applied during the GTLIKE 
routine, the 7-ray emission that is associated with Pictor 
A is located at (aj2ooo, <5j2ooo = 80.17°, —45.61°), with a 
95% error radius of 0.08°. The GTFINDSRC routine was 
also able to identify and localise the emission from both 
BZQ J0515— 4556 and the other TS 'hot-spot' as two 7-ray 
sources separate from Pictor A. 

Once satisfied that all point sources within an 8° ROI 
have been accounted for, we applied GTLIKE to the three- 
year data set. With the combined diffuse and 8° ROI point 
source model, we obtain the following best-fit power-law 
function for Pictor A: 

^ = (4.23 ± 0.44) x^^r-" " 

photons cm -2 s _1 MeV -1 (1) 

which equates to an integrated flux, in the 0.2 < E 1 < 300 
GeV energy range, of 

Fe>o.2 Gov = (5.8 ± 0.7) x 10" 9 photons cm" 2 s" 1 (2) 

only taking statistical errors into account. Primarily gov- 
erned by the uncertainty in the effective area, the systematic 

6 The test statistic, TS, is defined as twice the difference between 
the log-likelihood of two different models, 2[logL — logLo], where 
L and Lq are define d as the likelihood wh en the source is included 
or not respectively llMattox et al. Nl996tn . 
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Figure 1. Light curve of the 0.2 < E y < 300 GeV 7-ray flux 
from Pictor A during the first three years of Fermi-LAT observa- 
tions. The data is binned in yearly intervals, with the fluxes being 
plotted when the TS value is greater than 10 (~ 3<r), otherwise 
the 95% confidence upper limit is given. The bold horizontal line 
is the integrated flux given in Equation 2, derived from the three- 
year power-law fit. The difference in flux levels between the third 
annual bin and the three-year averaged flux, as well as between 
the second and third year, suggests variability on timescales of a 
year or less, though it should be stressed that further observations 
are required to confirm this. 

uncertainty of the integrated flux is energy dependent and 
is currently estimated as 10% at 100 MeV, do wn to 5% at 
560 Me V and back to 10% for 10 GeV photons (lAbdo et all 
<2011bh ). 

From the power-law fit, the predicted 7-ray count was 
337, with a test statistic of TS = 33.4, corresponding to a 
~ 5.8er detection. Compared to other 7-ray BLRGs, 3C 120 
and 3C 111, Pictor A's 7-ray flux is a factor of six lower, 
and the spectral index is comparable, with a spectral index 
of -2.93 ± 0.03 for Pictor A compared to -3.0 ± 0.3 for 3C 
120 and -2.7±0.2 for 3C 111. On the other hand, the 7-ray 
luminosity for Pictor A is comparable to both 3C 120 and 
3C 111, \og(L~,) — 43.1 ergs s - 1 compared to 43.6 e rgs s" 1 
and 43.9 ergs s _1 respectively (|Kataoka et al.1 |201ll 'l'). 

To investigate the possibility of temporal variability in 
Pictor A's 7-ray flux, the three- year data set was binned into 
one year periods and analysed separately with GTLIKE. 
As with the spectral analysis, only 0.2 < E 1 < 300 GeV 
events were considered, with the appropriate quality cuts 
applied, and the same combined diffuse and 8° ROI point 
source model was used. Only the fluxes of temporal bins 
with TS > 10 (~ 3<j) were plotted, with the 95% confidence 
upper limit being calculated for the remaining temporal bin. 
The resultant light curve can be seen in Figure [TJ with the 
bold horizontal line indicating the three- year averaged 0.2 < 
< 300 GeV flux level. 

A closer inspection of the individual photon energies 
reveals that the highest photon energy detected from Pictor 
A during this three year data set has an energy of 109.8 
GeV. Additionally, photons with energies of 99.7 and 91.3 
GeV were also detected, both within 0.8° of Pictor A. The 
spatial distribution of these highest energy photon events 
was compared to the sky map to confirm that none were 



spatially coincident with either BZQ J0515— 4556 or the 13.8 
TS pe ak. Taking the approach outlined in iNeronov et al.l 
(|2010r ) . the chance probability of the second highest photon 
being clustered within 1° of the highest photon is ~ 0.107. 
That is to say, the chance probability that these two photon 
events are not clustered background events is » 89.3%. This 
value is ~ 2a and therefore the observed 'clustering' of > 90 
GeV photons around Pictor A is not statistically significant. 
It is worth noting that a GTLIKE likelihood analysis for 
all very high energy (VHE) 7-ray events, 90 < £ 7 < 300 
GeV, from Pictor A, utilising the same combined diffuse and 
8° ROI point source model, finds a comparable statistical 
significance to Pictor A being a VHE 7-ray source, with a 
TS value of 4.99 or ~ 2.2a. 



4 SPECTRAL ENERGY DISTRIBUTION 
MODELLING 

We now turn to the question of whether the observed 7- 
ray emission is associated with the hot-spots in the terminal 
shocks of the radio lobes or with the relativistic jet. Fermi- 
LAT has resolved 7-ray emission from the radio lo b es of t he 
low-power radio galaxy, Centaurus A (|Abdo et al.l (|2010aT 0. 
However, beyond Centaurus A, further evidence for 7-ray 
emission from AGN radio lobes is scarce. 

As mentioned in §1, leptonic models which predict 7- 
ray emission from the radio lobe hot-spots attribute the ob- 
served radio-UV emission to synchrotron radiation, primar- 
ily due to the observation of pola rised radiation from the 
hot-spots at these wavelen gths fee. iRoeser fc Meisenheimerl 
l|l987h ; lHarris et al.1 (|l998l )V Given the synchrotron origin of 
this radiation, the shape of the SED at radio to UV wave- 
lengths places strong restrictions on the broken power-law 
energy distribution of the electron population within the 
hot-spot. The 7-ray emission, in such a model, is a result 
of the IC process acting on a percentage of the synchrotron 
photons. 

Due to the close proximity of Pictor A, radio lobe hot- 
spot features have been individually resolved at X-ray, op- 
tical, infrared and radio wavelengthfl These observations 
have discovered several bright knots within Pictor A's ra- 
dio lobe, with the Western hot-spot (WHS) being the most 
prominent and thus most extensively studied. SSC mod- 
elling of the WHS spectral energy distri b ution (SED) has 
already been performed by IZhang et al.l ([2009) , who pre- 
dicted a GeV 7-ray flux level on the order of ~ 10~ n ' 5 
ergs cm -2 s _1 . This prediction is somewhat higher than the 
three-year averaged flux that we have determined in the pre- 
vious section. 

Given there are a number of tuneable parameters avail- 
able in the SSC model we performed SSC modelling to 
determine the consistancy of the observed 7-ray flux be- 
ing produced through IC in the radio lobe hot-spots. The 
SED was constructed using historical radio, infra-red, op- 
tical and X-ray observations t hat have reported flux mea - 
surements from the WHS alone tMeiscnhcim er et al.l (|l997h ; 



7 The angular resolution of Fermi does not allow us to resolve 
Pictor A's radio lobes and nucleus regions separately. 
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Figure 2. SSC modelling of the SED of Pictor A's WHS 
from radio to 7-ray, plot t ed as uF v . The optical data is taken 



from [Mciscnhcimeretal 
taken from Tingav et al 



1 1997), the radio and infrared data is 
I 2008}), the X-ray data is taken from 
iMigliori et alj (120071) ? the VHE 7-ray upper limit, marked by an 
arrow, is taken from lAharonian et al.l ||200S| ) and the Fermi data 
point is report here. The solid line is the synchrotron spectrum 
and the dashed line is the inverse Compton spectrum. The SSC 
fit has been optimised to describe the radio to optical and the 
7-ray flux. 




Figure 3. SSC modelling of the SED of Pictor A's WHS 
from radio to 7-ray, plot t ed as vF v . The optical data is taken 

the radio and infrared data is 
the X-ray data is taken from 



from [Mciscnheimeretal 



taken from 



Tinga v et al 



[ 1997 ) 
2008) 



IMigliori et al. I (I2007T) . the VHE 7-ray upper limit, marked by an 
arrow, is taken from lAharonian et al.l ||2008| ) and the Fermi data 
point is report here. The solid line is the synchrotron spectrum 
and the dashed line is the inverse Compton spectrum. The SSC 
fit has been optimised to describe the radio to X-ray flux. 



iTineav et all (|200St l: IMigliori et all |2007h ). We also in- 
cluded the VHE 7-ray upper li mit set by the HESS collabo- 
ration l)Aharonian et all (2008)). The SSC model assumed a 
spherical emission region of radius R, moving with Doppler 
factor S and filled with a random magnetic field of strength 
B and an electron populat ion described by a broken power- 
law (for more details, see iKrawczvnski et all J2T11 1 1 i 1 . The 
radius of the emission region was fixed to the physical size 
of the WH S, R = 209 pc, a s obser ved with the Hubble Space 
Telescope (|Thomson et al.l (|l995l )) and the Doppler factor, 
S, was fixed to 1 thus removing any doopler boosting effects. 

In our first approach, the SSC fit to the WHS SED, 
shown in Figure [5] was originally optimised to account 
for the synchrotron radiation distribution, whilst accurately 
modelling our F ermi-LAT observations. The WHS SED was 
found to be best described by an electron distribution of 
Emin = 10 8 eV, E break = 4 x 10 11 eV and E max = 10 13 
eV, with broken power-law indices of ni = 2.4 and 712 = 4.4 
below and above the break energy respectively. To repro- 
duce the relative synchrotron and inverse Compton peak 
flux levels, a magnetic field of 4 x 10~ 5 Gauss and an elec- 
tron energy density of 2.5 x 10 -6 ergs cm -3 was required. 
The resultant SSC fit can be seen in Figure [2] with the solid 
line representing the synchrotron spectrum and the dashed 
line representing the IC spectrum. 

As can be seen in Figure[21 while the best fit SSC model 
adequately describes the radio-optical and Fermi fluxes, it 
over-predicts the inverse Compton flux at X-ray energies 
by an order of magnitude. It is worth noting at this point, 
that while smaller, this discrepancy is also present if the 
entire Western r adio l obe's X-ray flux is included in the SED 
(|Migliori et al.l (|2007l )). 

To address the over-prediction of the X-ray flux, both 
the magnetic field strength and electron population's en- 
ergy distribution were varied, optimising the fit for the syn- 



chrotron radiation distribution, whilst accurately modelling 
the X-ray observations. The best fit for the X-ray data 
was described an electron distribution of E m i n = 10* eV, 
E break = 2 x 10 11 eV and E max = 4x 10 13 eV, with bro- 
ken power-law indices of Tii = 2.4 and 712 = 4.4, a magnetic 
field of 1.5 x 10~ 4 Gauss and an electron energy density of 
2.5 x 10 -7 ergs cm -3 . It should be noted that a slight change 
in the break and maximum energies of the electron popula- 
tion were required so as to not over predict the flux of the 
high energy synchrotron emission. The resultant 'X-ray' fit 
can be seen in Figure [3] The fit parameters for Figure [2] and 
Figure [3] are shown in Table [TJ along with the fit parame- 
ters for Zhang et al.'s modelling of Pictor A's WHS and the 
Fermi-LAT Collaboration's modelling of Cen A's core and 
radio lobes. 



5 DISCUSSION 

In both Figure [2] and Figure [3] the radio to optical emis- 
sion can be seen to be accurately described by the process 
of synchrotron radiation. It is not, however, possible to si- 
multaneously describe the X-ray and 7-ray fluxes by the IC 
component of the one-zone SSC model. Therefore, our SED 
modeling provides support for the synchrotron origin of the 
low energy emission of the WHS of Pictor A. However, as- 
suming that the SSC model is the dominant mechanism for 
the observed X-ray flux from the WHS, then our SED mod- 
elling suggests that Pictor A's observed 7-ray emission is 
either produced through processes other than SSC, or that 
the 7-ray originates from either Pictor A's relativistic jet or 
central core region. 

Other possible mechanisms through which the 7-rays 
can be emitted include the external-Compton (EC) model or 
hadronic models invoking photon-proton or proton-proton 
interactions. Similiar to the SSC model, the EC model sim- 
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Table 1. Summary of fit parameters. The Cen A core fit paramete rs are taken from lAbdo et al. (2010a) while the 
Cen A radio lobe fit parameters are taken from lAbdo et al.l fcoioj l. It should be noted that the Cen A radio lobe 
fit considers 4 separate emission regions, 2 for each radio lobe; we quote the fit parameters for the total northern 
lobe, since the particle energetics for the two emission regio ns in the north lobe are the same. The Zhang et al. Pic 
A fit parameters are taken from the SSC fit parameters in lZhang et al. The energy limits of the electron 

population (E m ini Ebreak & E ma x) from our fits have been converted to Lorentz factors ("(mini 76reafc & 7mm) 
for comparison to the Cen A fits. 





Pic A (X-ray fit) Pic A (Fermi fit) 


Pic A (Zhang et al.) 


Cen A core 


Cen A !oije 


Magnetic field (G) 


1.5 x 10~ 4 


4 x 10~ 5 


4.4 X 10~ 5 


6.2 


8.9 x 10" 4 


Doppler factor (<5) 


1 


1 


1 


1 




volume (pc 3 ) 


2.14 x 10~ 7 


2.45 x 10~ 8 


4.1 x 10~ 7 


3.8 x 10~ 9 


3.6 x 10 15 




2.4 


2.4 


2.38 


1.8 


2.1 


r 2 


1.4 


4.4 


3.9 


4.3 


3.0 


7mm 


196 


196 




300 


1 


'y break 


7.8 x 10 5 


3.9 x 10 5 




800 


3.6 x 10 4 




2 x 10 7 


7.8 x 10 7 




1 x 10 8 


3.3 x 10 5 



ply has an external photon field, such as the Cosmic Mi- 
crowave Background (CMB), as the dominant photon field 
on which the inverse-Compton process acts. It is worth not- 
ing that Migliori et al. attributed t he X-ray emiss i on fro m 
the WHS to the EC of the CMB (jMidiori et all \200H )). 
however more recently Zhang et al. fo und the EC contrib u- 
tion from the CMB to be negligible (jZhang et all l|2009l )). 
While we have not considered EC during our SSC modelling 
of Pictor A's WHS, doing so would simply allow us to relax 
the energetic requirements of the electron population in the 
hot-spot. 

In depth hadronic interpretation of the Fermi-LAT ob- 
servations is outside the scope of this paper. Nonetheless, 
we do note that neutrino telescopes, such as ANTARES 
and IceCube, have the potential to constrain the hadronic 
contribution to Pictor A's observ ed 7- r ay flux (e.g. 

(1201 In ; [Abbasi et all 



iBecker et all (12007!)'): lAbbas et al 
|2012i ); I Adrian-Martinez et all (|2011 



))• 



However, it seems more likely that the reason why our 
SED modelling was unable to marry the WHS observa- 
tions with the observed 7-ray flux is that the 7-ray emis- 
sion is of jet origin. As mentioned in §1, 7-ray emission 
from BLRGs is traditionally interpreted as originating from 
a relativistic j et viewed at interme d iate a n gles from the line 
of si; ;h t (e.g. iGrandi fc Palumbol l|2007h ; ISambruna et al.l 
(2009); iKataoka et al l (|201ll )). Flux variability is a funda- 
mental property of astrophysical jets. Indeed, evidence for 
variability in Pictor A's jet has been observed with Chandra 
which has observed the X-ray flux from Pictor A's weste rn 
jet to vary on timescales of years ([Marshall et all (|2010h ). 
Evidence for variability in Pictor A's 7-ray flux can be 
seen in Figure [TJ which hints at 7-ray flux variability on 
timescales of a year or less. This evidence is primarily seen 
in the third year of observations, with a ~ 4a detection of 
Pictor A and a flux level that is 50% larger than the 3 year 
averagfl Taking into consideration statistical errors, the dif- 



8 It is also worth noting that the flux in the third year is approx- 
imately 30% larger than the 95% flux upper limit set in the first 
year of observations. Furthermore, the 3cr flux in the second year, 
is over five times smaller than the flux observed in the third year 
of observations. 



ference between the third year and three-year averaged flux 
levels is significant at the ~ 4a level. However, this is just 
short of the 5a level, and as such, further Fermi-LAT obser- 
vations are need to confirm this observational characteristic. 

At this point, we briefly turn our attention to other 
possible locations for the 7-ray emission, namely the core 
region of Pictor A within the immediate vicinity of the 
super- massive black hole (SMBH). Magnetospheric emis- 
sion models of AGN core regions have been successful in- 
voked to explain the 7-ray emission from other radio galax- 
ies such as M87 ( e.g.lRieger fc Aharoniaiil ((2008)). However, 
as highlighted in lBrown fc Adamsl (|201ll ). for these models 
Lie oc Mbh is expected, where Lie is the IC luminosity, 
which in this case is the 7-ray luminosity and Mbh is the 
SMBH mass. Given that Pictor A's SMBH mass is an order 
of magnitude larger then Cen A's, whilst the 7-ray flux is 
two orders of magnitude smaller, we find it unlikely that a 
magnetospheric emission model significantly contributes to 
the observed 7-ray emission and suggest that the observed 
7-ray emission from Pictor A is indeed of jet origin. 

Therefore, given that Chandra has resolved X-ray emis- 
sion from the WHS, the SSC fit of Figure [3] indicates that 
the observed 7-ray flux cannot be attributed to an SSC pro- 
cesses within the hot-spots of Pictor A's radio lobes. This 
conclusion, coupled with the evidence of variability observed 
in Figure [1] suggest that the 7-ray emission from Pictor A 
originates from within its jet, in agreement with the SED 
modelling of other 7-ray loud BLRGs. 



6 CONCLUSIONS 

We have reported the discovery of Pictor A as a 7-ray source 
with a significance of 5.8<r, utilising the first three years of 
observations with the Fermi-LAT detector. These observa- 
tions have found the E > 0.2 GeV spectrum to be well de- 
scribed by a power- law with V — —2.93 ±0.03, which is com- 
parable to the other Fermi-LAT detected BLRGs, 3C 120 
and 3C 111. Pictor A's three-year averaged 0.2 < F 7 < 300 
GeV flux is (5.8 ± 0.7) x 10" 9 photons cm" 2 s _1 , which 
equates to a 7-ray luminosity of log(L 7 )= 43.1 ergs s . 
Pictor A's 7-ray flux was thus found to be approximately 
six times smaller than 3C 120 and 3C 111, but due to its 
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close proximity, a luminosity that is comparable with other 
7-ray loud BLRGs. 

Pictor A's year-binned light curve suggested variability 
in the 7-ray flux, with the third year flux level being 50% 
larger than the three year average. However, further obser- 
vations are needed to confirm this observational property. 

SSC modelling of the WHS SED found that it was only 
possible to account for the 7-ray emission with an SSC model 
at the expense of greatly over-predicting the flux at X-ray 
energies. Given that the X-ray emission in the WHS has been 
resolved with the current suite of X-ray detectors, we suggest 
that the 7-ray emission from Pictor A originates from within 
its jet, in agreement with other BLRGs detected with Fermi. 
This suggestion agrees well with the evidence that the 7-ray 
flux is variable on timescales of a year or less. 

The highest energy photon detected during this 3 year 
period was 109.8 GeV, with several >90 GeV photons also 
being observed. A likelihood analysis of the 90 < F 7 < 
300 GeV 7-ray events does not find any significant VHE 
emission. While the clustering of the highest energy events 
does not appear to be significant, it raises the interesting 
possibility of VHE 7-ray emission from Pictor A. 
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